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Anisotropic electron-phonon interaction is shown to lead to the anisotropic polaron effect. The 
resulting anisotropy of the polaron band is an exponential function of the electron-phonon coupling 
and might be as big as 10 3 . This also makes anisotropy very sensitive to small changes of coupling 
and implies wide variations of anisotropy among compounds of similar structure. The isotope effect 
on mass anisotropy is predicted. Polaron masses are obtained by an exact Quantum Monte Carlo 
method. Implications for high-temperature superconductors are briefly discussed. 

PACS numbers: 63.20 Kr, 71.38+i 
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In an anisotropic crystal electron-phonon interaction 
(EPI) is anisotropic too. This is of minor importance 
in a metal with a high carrier density and weak EPI. 
The carrier mass is renormalized as to* = mo Q (l + 7 a A), 
where A is the electron-phonon coupling constant and 
7 Q ~ 1 a numerical coefficient corresponding to direction 
a. At A <C 1 mass renormalization is small and EPI has 
little effect on the band anisotropy which is governed by 
bare anisotropy of the rigid lattice. The situation is qual- 
itatively different in a semiconductor with a low carrier 
density and strong EPI A > 1. Under these conditions 
formation of small (lattice) polarons is expected |]-|. 
A general property of small polarons is an exponential 
renormalization of mass to* oc exp(7 a A/u>), where u> is 
the dimensionless phonon frequency GJ. Now the mass 
anisotropy rn^/mp oc exp((7 Q — r yp)X/Q) is an exponen- 
tial function of A, Q, and the form of EPI (through 7 a s). 
This simple relation hints several possible effects, (i) 
Anisotropy may be very large if A > 1 and Q < 1. (ii) 
Anisotropy is very sensitive to small changes of A caused, 
e.g., by doping (due to changes in screening). This 
also implies that different compounds that are similar in 
structure and have close A may nevertheless have very dif- 
ferent anisotropies. (iii) The dependence on phonon fre- 
quency is strong, hence the isotope effect on anisotropy. 
Thus anisotropic EPI may have profound effect on elec- 
tronic properties of doped semiconductors. 

The anisotropic case has been little touched in the po- 
laron literature. Kahn considered the large (continuum) 
Frohlich polaron with bare anisotropy and isotropic EPI 
and found that the polaron is less anisotropic than the 
bare carrier |j| . Recently Caprara and Del Prete found an 
enhancement of anisotropy by a factor «2 in the vicinity 
of a saddle point of the bare spectrum M. Both stud- 
ies were weak-coupling perturbational. In the strong- 
coupling regime the Lang-Firsov transformation |?J leads 
quite naturally to a strongly anisotropic polaron but we 
are not aware of any specific applications of this method. 
Also we are not aware of any previous numerical studies 



of the anisotropic case. Almost all the recent activity in 
the polaron field has been directed to the paradigmatic 
Holstein model || where EPI is local and the effect is 
absent. 

The recently developed Quantum Monte Carlo (QMC) 
algorithms ||,[l0| have opened the possibility to calculate 
the polaron effective mass exactly. The methods also 
allow studies of infinite lattices and arbitrary forms of 
electron-phonon interaction including anisotropic. In this 
paper we apply the method of Ref . jl(| , based on Feyn- 
man's integration of phonons |ll];[l^] , to two models with 
anisotropic EPI. The first one is a simple two-dimensional 
(2D) model where the enhancement of anisotropy is weak. 
The second model is a three-dimensional (3D) one with 
long-range EPI. Here we will find a very large enhance- 
ment and confirm the existence of the effects mentioned 
in the beginning. 

Let us begin with the simplest possible 2D model with 
anisotropic EPI (see Fig. |I|) . The lattice consists of two 
interpenetrating square sublattices n (crosses) and m 
(circles) . The sites n are fixed in their positions but sites 
m can vibrate along z-direction (up-down) , £ m being the 
internal coordinates. All m are uncoupled and have the 
same frequency uj and mass M. The carrier moves in 
the sublattice n by hopping between nearest neighbors 
with amplitudes i nn ' . The carrier at a site n attracts the 
two nearest sites m (above and below) with force n. The 
model Hamiltonian is 



H 
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where p m = —ih~d/d£ m . The hopping integral in x and z 
directions are t and t±, respectively. The model is char- 
acterised by the dimensionless frequency uj = Huj/t, the 
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FIG. 1. The two-dimensional lattice (positioned in the 
plane xz) of the model (Q). The carrier (bullet) hops be- 
tween n sites (crosses) but interacts with m sites (circles) 
which vibrate along z-direction. Two m sites are shifted from 
their equilibrium positions by some distance forming a small 
polaron. The three-dimensional lattice of the model (Q) is 
obtained by repeating the whole figure along y-direction per- 
pendicular to the sheet plane. 

bare anisotropy ratio rno z /mQ X = t/tj_, and the coupling 
constant A = Em /m(0)]/(2Mw 2 L>) = 2k 2 /{2Muj 2 D), 
where D is half of the bare bandwidth. If the carrier in- 
teracted only with one neighboring m site then Eq. (|^) 
would have described the 2D Holstein model. Polaron 
masses of the latter model for w=1.0 and t± = 0.5 1, ob- 
tained with QMC, are shown in Fig. ^ (main picture, 
squares). At A>1.5 both to* and to* increase expo- 
nentially as ocexp(4.2A) . However, the ratio to*/to* 
remains constant as a function of A and equal to the 
bare anisotropy 2 (see inset, squares). Thus the local 
interaction of the Holstein model does not enhance bare 
anisotropy. 

The behavior of the model (|J) is different. While the 
renormalization of both masses is similar to the Holstein 
case to* now grows faster than to* even for ij_ = t. The 
ratio to*/to* grows exponentially as oc exp(0.58 A) in the 
small-polaron regime A > 1.2 (see the inset in Fig. [|, 
circles). The physical reason for the difference is as fol- 
lows. In the small-polaron regime the carrier is mostly 
confined to a single lattice site n and both nearest m 
ions are shifted by some distance d from their equilib- 
rium positions (see Fig. |l|). Tunneling along axdirection 
occurs when two neighboring m ions on the right (or left) 
get shifted fluctuationally by d which creates a potential 
well on the neighboring n site, similar to the one the 
electron is already in. However, for the tunneling along 
z-direction to become possible only one m ion must be 
shifted by d but another, the one which is forming the 
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FIG. 2. Polaron masses in the model (|^) for t± = t (cir- 
cles) and in the 2D Holstein model for t±/t = 0.5 (squares). 
uj = 1.0 in both cases. Open (filled) symbols are m* (mj). 
Inset: anisotropy m*/m* for the model (Q) (circles) and Hol- 
stein model (squares). 

current potential well, must be shifted by the double dis- 
tance 2d. This process is less probable and requires more 
time to occur. That is why, although the lattice defor- 
mation extends in z-direction, it is more mobile in in- 
direction. 

The simple model considered illustrates the two generic 
properties of strong anisotropic electron-phonon inter- 
action: (i) it makes the spectrum anisotropic, and (ii) 
anisotropy is an exponential function of coupling. How- 
ever, the actual enhancement is modest in this model, 
to* /to* « 2 for A = 2.0 when the masses themselves are 
already huge, to* ~ 10 3 . In order to understand un- 
der which conditions anisotropy can be larger one has 
to consider how the form of electron-phonon interaction, 
i.e. the shape of the function / m (n) affects mass enhance- 
ment of the small polaron. Within the path-integral ap- 
proach the inverse polaron mass is the second moment of 
the end-to-end distribution of imaginary-time paths with 
open boundary conditions |l0[] . Each path includes in its 
statistical weight a phonon-induced factor e A where A 
is the action of the retarded self-interaction p| . Essen- 
tially, A depends on the forces / m (n) via combinations 

$(r 1 -r 2 ) = ^/ m (r 1 )/ m (r 2 ). (3) 

m 

Here r is the carrier position on a path, i.e., one of n. 
Normally, the function $(r) is the largest at r = and 
decays monotonically at large r. What matters for po- 
laron mass are the gradients V$(r) in the corresponding 
direction. The smaller the gradient the less the action 
changes upon the increase of ri — r2 , the more likely such 
a change will be accepted by the Monte Carlo process, 
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FIG. 3. Function «T 2 $(r) [Eq. (§)] for: (a) - the 2D model 
(|); (b) - the 3D model (0). 



the larger the second moment of the resulting end-to-end 
distribution. Thus, a steep (smooth) $(r) leads to a large 
(small) effective mass. 

For the Holstcin model (for which / m (n)=K(5 mn ) func- 
tion (||) is the lattice delta function 3>(ri— ra) = t 2 <5 ri r 2 
which is a steep function. Therefore in the small-polaron 
regime the weight of paths with non-zero ri — r2 is expo- 
nentially small leading to an exponentially small second 
moment. For the model (||) the x- and z- cross-sections 
of <&(r) are shown in Fig. ||(a). One can see that due 
to negative values at (x = 0, z = ±1) the gradient in z 
direction is bigger than in x direction. This is the reason 
for the enhanced to*. There are cases where V$(r) in 
one direction is much smaller than in the Holstein model 
or than in model (0). Consider again a lattice with the 
topology of Fig. nT Elementary geometrical considera- 
tions show that <&(!") is small when (i) the distance 
between chains n and m is large; (ii) interaction with 
distant m ions is included. When these conditions are 
fulfilled one should expect a much lighter cc-mass than in 
the Holstein case. At the same time V z <I>(r) is always big 
at the origin. Indeed, by definition $(0, 0) is positive but 
$(0, ±1) is negative (because / m (0, 0) and / m (0, 1) have 
different signs for m ions with m = (x, 1/2)). Therefore 
the z-mass is of the order of the Holstein mass and, con- 
sequently, a large anisotropy m* z /m* might follow. One 
should add, that if m ions were allowed to move along re- 
direction it would have increased m* and decreased m*. 
A predominant interaction with a single phonon polariza- 
tion is essential for a large anisotropy. Finally, it is not 
difficult to see that in 3D the difference between V 2 <I>(r) 
and V x Q(r) would be even bigger than in 2D. 

We now consider a 3D model where all the above 
factors are present. The lattice structure is obtained 




FIG. 4. Mass anisotropy of the model (|4|). Circles - 
Gj = 1.0, t± — 0.25 1; squares - u> = 0.5, t± — 0.25 1; tri- 
angles - Q = 0.5, t± — 0.1 1. Inset: masses m% y (filled circles) 
and m* for ui = 1.0, t± — 0.25 1. 



from the one depicted in Fig. g, by repeating it along y- 
direction which is perpendicular to the sheet plane. The 
distance between neighboring n and m planes is equal to 
the lattice constant in x and y directions (i.e., the lattice 
constant in z direction is twice the one in xy plane) . The 
EPI has the form 



/m(n) 



13/2 



(4) 



which is the z-projection of the Coulomb force with a 
screening radius R. This model is designed to describe 
the layered structure of high-T c superconductors, n sites 
(crosses) representing the copper-oxygen planes and m 
sites (circles) the apical oxygens [|l3| . Its ID (one n chain 
and one m chain) and 2D (one n plane and one m plane) 
versions were considered in p3] (see also |[o| ) . There was 
found a much reduced polaron mass, as compared with 
the Holstein model. Here we consider the full 3D prob- 
lem. Due to the low c-axis conductivity of the cuprates 
screening in z direction is poor and R is large. We use 
R = 10 lattice constants in x direction. For this interac- 
tion A = 2.93k 2 /(2Mui 2 D). Function $(r) for this model 
is shown in Fig. H(b). 

Effective masses of model (0J for ui = 1.0 and t± = 
0.25 1, obtained with the exact QMC, are shown in the 
inset of Fig [|. There is a striking difference in A- 



dependences of m* y and 



While 
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grows much 



slower than the Holstein mass, which is consistent with 
previous findings in one and two dimensions, to* grows 
similar to the Holstein case. For example, for A = 2.0 
to* = 420TOOZ while to* = IOtooz, for A = 2.5 m* = 
5600 TtiQ X while m* = 2lmQ X , etc. As a result, the mass 
anisotropy increases sharply with coupling reaching 2500 
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for A = 3.0 (see the main picture in Fig. ||, circles). Ba- 
sically, we are dealing with a new situation when two 
types of polaron are present at the same time. In z di- 
rection the carrier is the Holstein small polaron with a 
small size of lattice deformation (of the order of one lat- 
tice constant), localization of the carrier within this lo- 
cal deformation, and a very large effective mass. In xy 
plane the lattice deformation is larger than both the lat- 
tice constant and the localization area of the carrier, and 
the polaron is much more mobile (such a quasiparticle 
was named small Frohlich polaron in pj]). It is inter- 
esting to note fits m*/mo x — exp(1.26A + 0.88A 2 ) and 
m* z /m* y = exp(0.42A + 0.71A 2 ). The quadratic terms in 
the exponents are unusual. We were unable to establish if 
this is the true asymptotic behavior or the dependencies 
approach a pure exponential growth at larger A. 

We continue to study the model (f|) by changing its 
parameters. A general property of the small polaron 
is that its mass increases as the phonon frequency de- 
creases because the lattice deformation becomes less mo- 
bile, to* oc exp(7 ct A/[D). In our case j z > j xv . Therefore 
one should expect much stronger change of to* than of 
m* y with frequency. Figure^ shows m* z /m* y for Co = 0.5 
and tj_ = 0.25 t (squares). One can see that in this 
case the anisotropy is much bigger than when Co = 1.0. 
This implies the existence of the isotope effect on mass 
anisotropy in systems with strong long-range anisotropic 
electron-phonon interaction. Another model parameter 
is the bare anisotropy tj_/t. It has been found to be 
of minor importance. In short, it is a constant factor 
which is carried throughout the whole A interval. Fig- 
ure ^ shows m* z /ra* xy for Q = 0.5 and t± = 0.1 1 which 
corresponds to the bare anisotropy mo z / mo xy = 2.5 (tri- 
angles). The plot is very close to the one with t±_ = 0.25 1, 
the small difference resulting from the latter case being 
more adiabatic than t± = O.lt. 

Layered high-T c superconductors seem to be a good 
candidate for the effect described. Indeed, in the cuprates 
the carrier density is low, EPI is obviously anisotropic 
and strong, as revealed by isotope substitution [ ^5|Jl6| 
and neutron scattering flljj experiments. At the same 
time, a number of unusual c-axis properties fl8|] 
anomalously low interlayer hopping (large to*), anoma- 
lously large mass anisotropy, large variation of anisotropy 
among different cuprates — are reproduced in our model 
calculation. Although other effects, such as the pair- 
ing of carriers, have been shown to account for some c- 
axis anomalies of the cuprates the effect of strong 
anisotropic EPI seems to be worth further studying. 
For instance, the puzzling exponential increase of c-axis 
plasma frequency with doping and the exponential de- 
crease of anisotropy with doping observed in YBCO and 
other cuprates 1 18 can be explained easily by the present 
mechanism if one assumes that extra holes screen the 
long-range EPI and thereby reduce A ]l9| . 

In summary, it has been argued that strong anisotropic 



electron-phonon interaction in general enhances the 
anisotropy of the carrier band. The effect is strongest 
when (i) the interaction is predominantly with one 
phonon polarization only; (ii) the interaction is long- 
range; (iii) the dimensionality is three; (iv) the distance 
between planes is large. The mass anisotropy is an ex- 
ponential function of coupling which results in the large 
value of anisotropy and in wide variations among differ- 
ent compounds of similar structure. The bare anisotropy 
is irrelevant for the effect. With decreasing phonon fre- 
quency the mass anisotropy increases drastically. This 
amounts to isotope shift of the mass anisotropy. Possible 
candidates for the effect are layered high-T c supercon- 
ductors. 
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